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The multiply differential cross section of the ionization of hydrogen molecular ion by fast electron 
impact is calculated by a direct approach, which involves the reduction of the initial 6D Schrodinger 
equation to a 3D evolution problem followed by the modeling of the wave packet dynamics. This 
approach avoids the use of stationary Coulomb two-centre functions of the continuous spectrum of 
the ejected electron which demands cumbersome calculations. The results obtained, after verification 
of the procedure in the case atomic hydrogen, reveal interesting mechanisms in the case of small 
scattering angles. 

PACS number(s): 34.80.Dp 



I. INTRODUCTION 

New experimental methods, particularly, based on the 
multiple coincidence detection technique [0-^ stimulate 
the interest to fundamental theoretical studies of the dis- 
sociative ionization of diatomic molecules by electron im- 
pact. In this context the molecular hydrogen ion can be 
considered as the basic system in which the removal of the 
unique electron causes dissociation. Substantial theoret- 
ical analysis of the dissociative ionization of by fast 
electrons was recently carried out in [Q. As mentioned 
in the crucial point of calculating the cross-section 
of such processes is that no closed exact analytical wave 
functions of the continuum states exist. In Q the final- 
state wave function of the ejected electron was found by 
taking a product of two approximate functions that take 
into account the two scattering centers. To improve the 
calculation it seems straightforward to obtain these func- 
tions with the exact numerical solutions of the two-center 
continuum problem. However, this approach involves a 
cumbersome procedure of calculating multi-dimensional 
integrals of the functions presented numerically that re- 
quires huge computer facilities and may cause additional 
computational problems. It seems reasonable to search 
for direct computational approaches, in which the basis 
of exact two-center continuum wave functions is not in- 
volved. Note that the potential advantage of such meth- 
ods is that they could be generalized over a wider class of 
two-center systems starting from the molecular hydrogen 
ion as a test object. In the present paper we develop a 
direct approach to the ionization of hydrogen molecular 
ion by fast electrons that involves the reduction of the ini- 



tial 6D Schrodinger equation to a 3D evolution problem 
followed by modeling of the wave packet dynamics. 

Originally we intended to treat the incoming elec- 
tron classically, its trajectory being approximated by a 
straight line with the deflection neglected. The bound 
electron was to be treated quantum mechanically. Pre- 
liminary calculations at the impact parameter p = 10 
a.u. has shown, first, that the probability of the emission 
of the electron having the energy of 50 eV is extremely 
small, and, second, that the direction of the electron 
emission is orthogonal to that of the incoming electron 
motion, that contradicts the results of Q]. This means 
that the main contribution to the small-angle scattering 
comes from the central collisions with the bound electron 
in the region of its localization. Generally, the classically 
estimated deflection angle of 1° for scattered electron cor- 
responds to the impact parameter of the order of 1 a.u., 
so that the trajectory passes through the molecule and 
the classical treatment of the incoming electron is not 
valid. 

Here we develop and apply a direct approach to the 
calculation of the angular distribution of scattered and 
ejected electrons that involves the reduction of the ini- 
tial 6D Schrodinger equation to a 3D evolution problem 
followed by modeling of the wave packet dynamics. The 
approach does not make use of the basis of stationary 
Coulomb two-center functions of the continuous spec- 
trum for the ejected electron, whose proper choice is a 
crucial point of other model calculations. Our approach 
can be considered as the linearized version of the phase 
function method for the multi-dimensional scatter- 
ing problem. The evolution problem is solved using the 
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method based on the spHt-step technique with com- 
plex scahng, recently proposed by some of us and tested 
in paraxial optics ||]. In the present paper the method 
as a whole is also tested using the well known problem of 
electron scattering by hydrogen atom M. 



II. BASIC EQUATIONS 

We start from the 6D stationary Schrodinger equation 
which describes two electrons in the field of two fixed 
protons 



*(r,R) = £;*(r,R), (1) 




^ Y,y 

FIG. 1. Coordinate frame 

where r is the radius-vector of the electron initially 
bound in and finally ejected, R is the radius-vector 
of the impact electron, Hq — — ^\/-^+U{r) is Hamiltonian 
of ejected electron in the field of two protons, y(r,R) = 
?7(R) -I- C/i„t(r, R) is the interaction between the impact 
electron and molecular ion, U{r) = — 1/ri — l/r2 is the 
attractive potential between the ejected (scattered) elec- 
tron and the protons, ri = |r — rip|, r2 = |r — r2p|, 
Tip is the radius- vector of the i-th proton, C/i„t (r,R) = 
l/|r— R| is the repulsive potential of interaction between 
the electrons. The origin of the coordinate frame is cho- 
sen in the center of symmetry of the molecular ion with 
the Z axis directed along the momentum of the incident 
electron. 

For the scattering problem solved here the energy of 
the system may be presented as E = kf /2 + Eq, where 
—Eq is the ionization potential, ki is the momentum of 
the incident electron. Let us seek the solution of Eq.(j^) 
in the form ^{r,X,Y,Z) = V'lr, R±, Z) exp(ifcjZ). Un- 
der the condition that (fcg +k]_ — 2Eo)/kf « 1 one can 
neglect the second derivative of i/j with respect to Z. As a 
result we get the evolution-like equation for the envelope 
bmction ijj{r, Il±, Z) 



ik. 



dZ 



^o(r) - -V^^ --Eo + Vir,R) \ V(r, R^, (2) 



Neglecting the large- angle scattering one can write the 
initial condition for ip as 



■0(r, Ri,-oo) = ^o(r) 



(3) 



To solve the 5D Schrodinger evolution equation(^ we 
use Fourier transformation with respect to the variable 



Ri 



V'(r, Ri, = ^ J V'k^ (r, Z) exp(^kiRi)dRi. (4) 
Then Eq.(||) takes the form 

+ I ^k,k'^ (r, Z)V'k'^ (r, Z)dk'^, (5) 

where 

Vk^k'^(r,Z) = J exp(-z(ki-kl)Ri)I^(r,Ri,Z)dRi 

(6) 

is the Fourier transform of the interaction potential 
Vir,R^,Z). 

Further simplification of the problem is possible if the 
amplitude of the incident wave is much greater than that 
of the scattered wave. In this case one can put 



7/^kJr,Z) =^(ki)V^o(r) 



(7) 



in the integral term of Eq.(|5|). As a result we get the 
inhomogeneous equation 



dZ 



+ 
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-So Uk,(r,Z) 



(8) 



where Vk^ (r, Z) = Vk^o(r, Z), with the initial condition 
-0ki(r, -oo) = 0. 

To calculate the integral with respect to transverse 
variables in the expression for Vk^ (r, Z) it is easier to 
start from the known integral 



J exp(-ikR)-^dR : 
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Carrying out the inverse Fourier transformation 

dkz TT 



/ 

J — -( 



one gets 



exp{ikzZ) 



-k±\Z\ 



ki 



(9) 



(10) 
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^ — k±\Z-dz\-i'k-±d± _^ g— fej_ | Z+dz | +ikj_dj 



(11) 
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Here k± = ki sin 9s is the transverse momentum compo- 
nent of the scattered electron, Og is the scattering angle, 
±d are the positions of the nuclei with respect to the 
center of symmetry. Note that terms in square brackets 
determine the elastic scattering of the incident electron 
by the nuclei. 

Due to the exponential decrease of the source term 
with \Z\ the integration may be actually carried out 
within a certain finite interval {—Zmax, Zmax)- Hence 
the zero initial condition should be imposed at the point 

_ 7 

^max ■ 

Note that the approximation is actually equivalent 
to the first Born approximation^. Multiply Eq.(^) by 
the complex conjugate function of the continuous spec- 
trum of Hq and integrate over all r. Then 



ik. 



dCk^(ke,Z) 



dZ 



(27r) 



^*(ke,r)Fk^(r,Z)^o(r)dr, 



where Ck^ (kg, 2') — J V^*(ke, r)?/'kj_ (i": Z)dr is the proba- 
bility density amplitude for the transition of the initially 
bound electron into the state with the momentum kg. 
Let us substitute 

Ck^(ke,Z) = (7k^(ke,Z)exp(ifczZ), 

where fc^ is the increment of the longitudinal component 
of the momentum of the impact electron determined by 
the relation 



1 /fc2 ^.2 



(13) 



This relation is actually equivalent to the energy conser- 
vation law written neglecting the terms of the order of 
fc|. The substitution yields 



ik. 



1 



dZ 

exp{~ikzZ) / r(ke,r)Fk^(r,Z)V'o(r)dr, (14) 



(2^)2 
and 



Ck^(ke,CX)) = 

-^^(e^'^=^^(ke,r)|T/(r,R)|e^'''^Vo(r)), (15) 

where ks = k^ — K is the momentum of the scattered 
electron, K = {^kx, —ky, —kz) is the momentum trans- 
fer. 

Provided that the ejected electron has the momentum 
ke, the asymptotic form of the solution of Eq. (|l|) for the 
wave function of the scattered electron when i? — > oo is 

= eM^k^Z) + £^^P(!^/k^((?„0,). (16) 



The scattering differential cross-section(DCS) can be 
then expressed as 



(17) 



On the other hand, the asymptotic form of the wave 
function resulting from the solution of Eq. (^) under the 
condition Z oo can be presented as 

*^^(R) exp(ifc,Z)-f 

exp{ikiZ) / C-ki^Cke, oo) exp{i'k±IL± + ikzZ)d'k±. (18) 



Making use of the fact that the integrand has a stationary 
point we finally get 



iki Z 



1 



(-27^^fcOC'kO (ke,cx)), (19) 



where k']^ — fc^ sin 0s (cos 0s, sin 0^), R± = RsmOg, Z = 
R cos 9s . The expression (p9|) agrees with ( p^ ) within the 
accuracy of the order of 92ii we set 

fkA(^s,<f>s) = -2TrikiCy,o (ke,oo) 



2tt 



(e*'^=^^(ke , r)\V{v, R) le^'^'^Vo (r)) . (20) 



The latter expression is similar to the formula for 
/ke(^s,0s) derived in |^ using the first Born approxi- 
mation. 



III. CALCULATION OF THE ANGULAR 
DISTRIBUTION 

The asymptotic expression of the radial part of the 
wave function corresponding to the continuous spectrum 
of Hq can be written as 



1 



exp{-iEt + i / v{r')dr') (21) 



where t = Zjki is the evolution variable, v(r) — 

^2(£;-C/-(r)), £; = ^ + 4 - i^o, C/"^(r) = -Z' jr, 
Z' = 2 is charge of two protons. In the asymptotic limit 
one can take only the radial component of the momen- 
tum of the ejected electron into account. Then, accord- 
ing to pO| ], the expression for calculating the amplitude 
A{k, 9, ^takes the form 



Ak^(fce, 9e, 4>e) 



1 



where 



dt'j{4'kAr,0,^,t'),rE{r,t')) 



(22) 
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dr 



5* 

dr 



(23) 



is the flux introduced in pO| ], to = —Zmax/ki and ti >> 
Zmax/ki- The approximate relation ( p2[ ) becomes exact 
when ti — > +CXD and simultaneously r„iaj; — > +00. 

The amplitudes defined by (03) are related with the 
coefficients introduced in Eq.(|15D by 



\A]^^{ke,9e,4>e)\ = fcg Ckj^ (kg , OO 



(24) 



Using (|l7|),(|2ClD and (|2j) we get the final expression for 
the differential cross-section 



(25) 



In the region where r > r^ax we made use of the com- 
plex scaling technique pT| ] to suppress the non-physical 
reflection from the grid boundary. 



IV. NUMERICAL SCHEME 

The inhomogeneous Schrodinger equation to be solved 
can be written as 



a^'(r,t) - . . . . . . 

z^^=i7o(r)vI/(r,i)-fF(r,i), 



(26) 



The solution of Eq. (^6|) to within the second-order terms 
in At can be expressed as the following sequence of equa- 
tions: 



iAt , , 
^it)-—F{r,ty, 



iAt 



(27) 
(28) 



^'{t + At) = ~ —F{r,t + At). (29) 

The key step of the procedure is Eq. (|8|) which de- 
fines nothing but Cranck-Nicholson scheme. To solve this 
equation we make use of the partial coordinate splitting 
(PCS). A finite-difference scheme is applied for the radial 
variable r and the polar angle 9. Fast Fourier transform 
(FFT) is used for the azimuthal angle <j>. 

In the spherical coordinate system, the z-axis of which 
is directed along the symmetry axis of the molecule (and 
not along the velocity of the impact electron) and after 
the substitution ^ = ip/r and the Fourier transformation 



of Eq.(|28|) the terms Hq"^ 
into 

Ho{r,iri,m)'ip{r,r],m) 



0,1 



entering this equation turn 
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1 



5 n 2n 5 



drj 1 — if 



ip{r, rj, m) 
(30) 



+U{r, r])ij{r,r],m), 
where rj — cosO, m is the asimuthal quantum number 



and 



Finite-difference approximation 

yields M sets of linear equations, 

L and N 



D'' - of the differential opera- 



tors entering Eq.(C 

each set being of the order L x N, where M, 
are the numbers of grid points in 0, rj and r, respectively. 
Direct solution of each set of equations requires NL'^ op- 
erations at each step in t. The FFT that should be per- 
formed twice, first, when proceeding from ( ^7|) to (|2^), 
and second, from ( |2^ ) to (p9|), requires NLM log2 M ex- 
tra operations. 

To reduce the number of operation we propose a 
double-cycle split-step scheme. In case when Ho can be 
presented as a sum Hq = Hi + H2 , this scheme can be 
formulated as follows 





■01 




{1+ 


At ^ , 

1—Hi)'ljj2 


At ^ , 
- {I-i—Hi)^i 


(1+ 


i^H2)ljj3 


= (/-J^i?2)02 


(1+ 


i^H2)->jM 


= (/-Z^i?2)03 


{1+ 


i^Hi)->Jj5 


= {I-t^H^)^^ 




ip{t + At) 


= "05, 



which to within the second-order terms in At corresponds 
to the initial Cranck-Nicholson scheme, / and Hi, 2 is 
square matrixes, {I)i2iij2ji = 5i2ii5j2]i- 

Now the problem is how to split the Hamiltonian Hq 
into two parts. Formal separation of radial and angular 
parts leads to difficulties associated with the singularity 
of the angular part. Due to this singularity the scheme 
appears to be conditionally stable with severe limitations 
imposed on the step At. Practically this version of the 
splitting scheme is applicable only if the grid in r is rough 
enough. 

To remove this limitation we propose a partial coor- 
dinate splitting scheme. Its principal idea is that in the 
vicinity of r = it is preferable not to split off the angu- 
lar part at all. To implement this idea we introduce the 
r-dependcnt weight function p{r) which is supposed to 
diminish in the vicinity of r = and define the discrete 
approximation of the operators Hi 2 in the following way 



+p{n^) 

(-^2 )'i2ni2ii 



1 

''''' 5^^^,+U2{n,,Tjj,) 



2 rl 



1 D"'™ 



(31) 



(32) 



here Uas{'r)^U2{r, rf) = ?7(r, rf). It is reasonable to choose 
p(r) as a cubic polynomial 
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p(r) 



+ 1, ra < r < Ta + ap] 
r >ra + ap] 



where is the radius of the vicinity of r = where the 
sphtting is absent, Up is the width of the area of partial 
spHtting. Such a polynomial satisfies the condition of 
smooth connection at the boundaries that separate the 
region of partial splitting from the regions of full split- 
ting, on one hand, and of no splitting at all, on another 
hand. 



V. NUMERICAL CALCULATIONS AND 
RESULTS 

The method was tested using the well-studied exam- 
ple of the impact ionization of atomic hydrogen. We 
compared our results with those given by the well-known 



expression obtained in the first Born approximation [ [L2[ . 
Good agreement was demonstrated in the energy interval 
of interest from 1 to 3 a.u., E^, being the energy of 
the ejected electron. 




FIG. 2. The multi-fold differential cross section (MDCS) of 
the ionization of versus the ejection angle and ejection 
energy _Ee for Qd ~ 135° 




45 90 135 180 225 270 315 360 

(deg.) 



0.04 



(a) 




135 180 225 270 315 360 

(deg.) 



(c) 




135 180 225 270 315 360 



FIG. 3. The multi-fold differential cross section (MDCS) of the ionization of H2 versus the ejection angle 9^ for different 
angles Od: a)6ld = 139.2° that corresponds to d || K; h)ed = 49.2° that corresponds to d _L K ; c)9d = 0°; d) dd = 90°. The 
energy of the ejected electron Ee — 1.85 a.u. =50. 3 eV. 
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Our numerical studies concerning the molecular hydro- 
gen ion focused on the variation of the multi-fold differ- 
ential cross section (MDCS) concerning a coincidence de- 
tection of the two emerging electrons and one of the pro- 
tons with the ejection angle 9^ at different orientations of 
the molecular axis, provided that the scattering angle is 
small. The examples of our results illustrated by Figs.^j^ 
are obtained under the following conditions: the momen- 
tum of the impact electron ki =12.13 a.u. ~ 2000 
eV); the angle of scattering 9s = 1°. The impact and 
ejected electron trajectories and the molecular axis are 
supposed to lie in one plane. The latter restriction is 
not imposed by the method as such, it is just an ex- 
ample. Generally, one gets full information about the 
ejected electron, i.e., the dependence of MDCS from Ee, 
9e and 0e, after each run of the code at given values of 
the impact energy, scattering angle and molecular axis 
orientation. In Fig.H demonstrates the energy-angle dis- 
tribution, extracted from the data getting in result of one 
run of the code. In the planar geometry the orientation 
of the molecular ion is determined by a single angle 0^ be- 
tween the impact direction and the internuclear axis. We 
remind that the momentum transfer vector was defined 
above as K = — kg. In Figs. || we present the par- 
ticular cases of the dependence of MDCS upon 9e when 
internuclear axis is a)parallel to the momentum transfer; 
b)perpendicular to the momentum transfer; c)parallel to 
the impact electron direction k^; d)perpendicular to the 
impact electron direction k^. As it could be expected 
basing on the elementary symmetry considerations, the 
first two plots are symmetric with respect to the direction 
of the momentum transfer that corresponds to the angle 
9e = 319.2°. Since this symmetry is not assumed a pri- 
ori in the procedure, this may be considered as one more 
evidence in favour of the validity of the results demon- 
strated. 




FIG. 4. The multi-fold differential cross section (MDCS) of 
the ionization of H2 versus the ejection angle 6e and molec- 
ular angle Od- The energy of the ejected electron Ee = 1.85 
a.u.=50.3 eV. 

The recoil momentum Qrecoi/ = K — kg transmitted 
to the target has its minimum for ke parallel to K. In 
this case all the momentum is transferred to the ejected 



electron and the probability of the ionization is maximal. 
This is confirmed around 9e = 319.2° on figures ||(a) and 
^(b) where the inter-nuclear axis is respectively perpen- 
dicular and parallel to K. So this is a good verification for 
our calculation. Now, for the situation where kg is anti- 
parallel to K, the recoil momentum Q recoil is maximal 
and the probability of the ionization is maximal. This 
is also visible for 9e — 139.2°. Now for the directions of 
the internuclear axis other than 9d = 139.2° (where d is 
parallel to K) or Od = 49.2° (where d is perpendicular to 
K) the target does not respect the above analysis. This 
is due to the fact that the diatomic target behaves as an 
atomic target only for these two angles. The other situ- 
ations present interference patterns the minima of which 
move when 9d changes. 

Fig.^ shows MDCS versus the ejection angle 9e and in- 
ternuclear angle 9d- As one can see, this dependence has 
rather a complex behaviour. 



0.08 r 




15 30 45 60 75 90 105 120 135 150 165 180 



(deg.) 

FIG. 5. The multi-fold differential cross section of the ion- 
ization of H2 as a function of the angle 9d between the im- 
pact direction and the internuclear axis for fixed ejection angle 
Be — 319.2°. The energy of the ejected electron is i5e = 1.85 
a.u. 

To confirm the above dependence we show in Fig.|| 
a section of Fig.^ for fixed ejection angle 9e — 319.2° 
which corresponds to the case when the ejected electron 
direction is parallel to the momentum transfer vector. It 
presents a variation of the MDCS with respect to the di- 
rection of the inter- nuclear axis. It can be clearly seen 
that the maximal value of MDCS is achieved when the in- 
ternuclear axis is perpendicular to the momentum trans- 
fer direction that correspond to 9d = 49.2°. This result 
agrees with the hypothesis formulated in |0[ . 



VI. CONCLUSION 

We have developed a procedure which determines the 
multiply differential cross section of the (e,2e) ionization 
of hydrogen molecular ion by fast electron impact, us- 
ing a direct approach which reduces the problem to a 
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3D evolution problem solved numerically. Our method 
avoids the cumbersome stationary perturbativc calcula- 
tions, and opens the way for near future applications to 
the (e,2e) ionization of more complex atomic and molec- 
ular targets. 
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